Dynamic transitions between metastable states in a superconducting ring 
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Applying the time-dependent Ginzburg-Landau equations, transitions between metastable states 
of a superconducting ring are investigated in the presence of an external magnetic field. It is shown 
that if the ring exhibits several metastable states at a particular magnetic field, the transition from 
one metastable state to another one is governed by both the relaxation time of the absolute value of 
the order parameter T|^| and the relaxation time of the phase of the order parameter r^. We found 
that the larger the ratio t\^\/t^ the closer the final state will be to the absolute minimum of the 
free energy, i.e. the thermodynamic equilibrium. The transition to the final state occurs through a 
subsequent set of single phase slips at a particular point along the ring. 
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I. INTRODUCTION 

The fast development of experimental techniques 
makes possible the study of physical properties of samples 
with sizes of about several nanometers. The interest in 
such small objects is due to the appearance of new effects 
if the size of the system is comparable to some charac- 
teristic length. In case of superconductors it means that 
the sample should have one size at least of order the co- 
herence length ^ (such superconductors are called 'meso- 
scopic' superconductors). For example, only for hollow 
cylinders or rings of radius ~ ^ it is possible to observe 
the famous Little-Parks oscillationstl. 

The majority of the research in this area has been 
limited to the study of static or quasi-static properties 
of mesoscopic superconducting rings, disks and other 
geometriesH. It is surprising that the dynamics in such 
systems was practically not studied. However, it knows 
out that the dynamics is very important for systems 
which exhibit a series of metastable states and which may 
be brought far from thermal equilibrium. For such sys- 
tems the fundamental problem is to determine the final 
state to which the system will transit too (see for example 
Ref.l). 

In the present work, we present for the first time a de- 
tailed study of the dynamic transitions between different 
states in a mesoscopic superconducting ring. This sys- 
tem is a typical example of the above mentioned systems 
where a set of metastable states exist. We will show 
that the final state depends crucially on the ratio be- 
tween the relaxation time of the absolute value of the 
order parameter and the time of change of the phase of 
the order parameter. Our theoretical results exjolain the 
recent magnetization results of Pedersen et alu in thin 
and narrow superconducting Al rings. 

To solve this problem we present a numerical study 
of the time-dependent Ginzburg-Landau (TDGL) equa- 
tions. Therefore, our results will also be applicable to 
other systems (for example, liquid helium) where the dy- 
namics are described by these or similar equations. In 



our approach we neglect the self-field of the ring (which 
is valid if the width and thickness of the ring are less than 
A - the London penetration length) and hence the distri- 
bution of the magnetic field and the vector potential are 
known functions. The time-dependent GL equations in 
this case are 



(V-iA)2V'-f7/;(l-|^n+x,(la) 



A(p = div (Im(V/(V - iA)?/;)) 



(lb) 



where ■0 



o«0 



is the order parameter, the vector po- 



tential A is scaled in units $0/(271"^) (where $0 is the 
quantum of magnetic flux), and the coordinates are in 
units of the coherence length ^(T). In these units the 
magnetic field is scaled by Hc2 and the current density, 
J> by jo = c$o/87r^A^^. Time is scaled in units of the 
Ginzburg-Landau relaxation time tgl = 47rcr„A^/c^, the 
electrostatic potential, (p, is in units of c<I>o/87r^^A(T„ ((t„ 
is the normal-state conductivity), and u is a relaxation 
constant. In our numerical calculations we used the new 
variable U — exp{—i J Adr) which guarantees gauge- 
invariance of the vector potential on the grid. We also 
introduced small white noise x in our system, the size of 
which is much smaller than the barrier height between 
the metastable states. 

We consider u as an adjustable parameter which is a 
measure of the different relaxation times (for example, 
the relaxation time of the absolute value of the order 
parameter) in the superconductor. This approach is mo- 
tivated by the following observations. From an analysis 
of the general microscopic equations, which are based 
on the BCS model, the relaxation constant u was de- 
termined in two limbing cases: u — 12 for dirty gap- 
^nd u = 5.79 for superconductors 
However, this microscopic the- 
ory is built on several assumptions - for example that the 
electron-electron interaction is of the BCS form, which 
influences the exact value of u. On the other hand, the 
stationary and time-dependent Ginzburg-Landau equa- 
tions, Eqs. (la,b), are in some sense more general and 
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their general form does not depend on the specific micro- 
scopic model (but the value of the parameters, of course, 
are determined by the microscopic theories)l3. 

The paper is organized as follows. In Sect. II we 
present our numerical results for the solution of Eqs. 
(la,b) for strictly one-dimensional ring. In Sect. Ill we 
explain the obtained results in terms of different time 
scales of the mesoscopic superconductor. The effect of 
the finite width of the ring is studied in Sect. IV where 
we also discuss the influence of non-zero temperature and 
finite A. 



II. SUPERCONDUCTING RING OUT OF 
EQUILIBRIUM 



(for i = we have = 2pc/R). 
and we can write Eq. (2) as 



In this case p 
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2B^' 



(3) 



For this value of the flux the thermodynamical equilib- 
rium state becomes L^q = Int(<I>c/$o). For example, if 
i? = 10 we find L^q = 6 (see Fig. 1). The fundamental 
question we want to answer is: what will he the actual 
value of the vorticity of the final state! Will it be the 
thermodynamic equilibrium state or a metastable one? 
This answer will be obtained from a numerical solution 
of the time-dependent Ginzburg-Landau equations. 



If the width of the ring is much smaller than ^ and 
X± = / d {d is the film thickness) it is possible to con- 
sider the ring, of radius i?, as a one-dimensional object. 
At first we limit ourselves to the solution of the strictly 
one-dimensional Eqs. (la,b) which model the dynamics 
of a ring with small width. In our approach the vector 
potential is equal to A ~ HR/2, where H is the applied 
magnetic field. 
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FIG. 1: Dependence of the free energy of the one-dimensional 
ring (with R = 10^) on the applied magnetic field for different 
vorticity L. 



As shown in Refs.t§ll3 the transition of the supercon- 
ducting ring from a state with vorticity L = f\74'ds/2TT 
(which in general can be a metastable state) to a state 
with a different vorticity occurs when the absolute value 
of the gauge-invariant momentum p = V0 — A reaches 
the critical value 
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Let us, for simplicity, consider that the magnetic field 
is increased from zero (with initial vorticity of the ring 
L = 0) to the critical He where p becomes equal to Pc 




GL 



FIG. 2: Dependence of the order parameter (a), the super- 
conducting js (b) and the normal j„ (c) current density in 
the point where the minimal (thick curves) and the maximal 
(thin curves) values of the order parameter is reached as func- 
tion of time for u = 3, R — 10 (dotted curves) and R = 15 
(solid curves). For both rings we took At = 100, the time 
interval over which the magnetic field is increased to its value 
He (inset of Fig. 2(a)). 

For our numerical calculations we considered two dif- 
ferent rings with radius R = 10 and R = 15. These 
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values were chosen such that the final state exhibits sev- 
eral metastablc states and, in principle, transitions can 
occur with jumps in L which are larger than unity. For 
R = 10, AL may attain values between 1 and 6 and for 
R = 15 it may range from 1 to 9. Smaller rings have 
a smaller range of possible AL values. We found that 
larger rings than E = 15 did not lead to new effects 
and are therefore not considered here. In our numeri- 
cal simulations we increased H gradually from zero to a 
value He + AiJ (we took AH = 0.036ifc for i? = 10 and 
Ai? = 0.012i/c for R = 15) during a time interval At, 
after which the magnetic field was kept constant. The 
magnetic field A_ff and time Af range was chosen suffi- 
ciently large in order to speed up the initial time for the 
nucleation of the phase slip process, but still sufficiently 
small in order to model real experimental situations in 
which H is increased during a time much larger than 
u ■ TGL (for example for Al tgl{T = 0) ~ 10""s). A 
change of AH and At, within realistic boundaries, did 
not have an influence on our flnal results. 

From our detailed numerical analysis, we found that 
the vorticity L of the final state of the ring depends on 
the value of u and is not necessary equal to L^q. The 
larger u the larger the vorticity after the transition. The 
final state is reached in the following manner. When the 
magnetic field increases the order parameter decreases 
(see inset of Fig. 2(a)). First, in a single point of the 
ring a local suppression, in comparison with other points, 
of the order parameter occurs which deepens gradually 
with time during the initial part of the development of 
the instability. This time sc;ale is taken as the time in 
which the value of the order parameter decreases from 1 
to 0.4 in its minimal point. When the order parameter 
reaches the value ^ 0.4 in its minimal point the process 
speeds up considerably and the order parameter starts 
to oscillate in time at this point along the ring (see Fig. 
2(a)). At the same time also oscillatory behavior of the 
superconducting js = Re(V'*(— iV — A)V-') (Fig. 2(b)) 
and the normal j„ = —Vip (Fig. 2(c)) current density is 
found at the point where the minimal value of the order 
parameter is found. At other places in the ring such 
oscillatory behavior is strongly damped (sec e.g. the thin 
curves in Fig. 2). After some time the system evolves to 
a new stable state, which in the situation of Fig. 2 is 
i = 3 for i? = 10 and L = 5 for i? = 15. 

In Fig. 3 the dependence on u of the dynamics of 
\tjj\ in the minimal point is illustrated. Note that with 
increasing value of u the number of phase slips (or equiv- 
alently the number of oscillations of the order parame- 
ter at its minimal point) increases and hence AL also 
increases. Another important property is that the am- 
plitude of those oscillations decreases with increasing u. 



III. TIME SCALES 

The above results lead us to conclude that in a su- 
peconducting ring there are two characteristics time 
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FIG. 3: The dynamics of the order parameter in its mini- 
mal point for different values of the parameter u and for two 
values of the radius R = 10^ (dotted curves) and R = 15^ 
(solid curves). The zero of time is taken at the moment when 
the value of the order parameter becomes equal to 0.4 in its 
minimal point. Nonzero values of near the moment of the 
phase slip is connected with the finite coordinate and time 
step used in the numerical calculations. 



scales. First, there is the relaxation time of the abso- 
lute value of the order parameter r|^|. The second time 
scale is determined by the time between the phase slips 
(PS), i.e. the period of oscillation of the value of the or- 
der parameter in its minimal point (see Fig. 3). Below 
we show that the latter time is directly proportional to 
the time of change of the phase of the order parameter 
T0 and is connected to the relax;ation time of the charge 
imbalance in the system. 

In Fig. 4 the distribution of the absolute value of the 
order parameter tp and the gauge-invariant momentum p 
are shown near the point where the first phase-slip occurs 
for a ring with radius i? = 10 and u = 3 at different times: 
just before and after the first PS which occurs at t ~ 
2-7tgl. Before the moment of the phase slip the order 
parameter decreases while after the PS it increases. In 
order to understand this different behavior let us rewrite 
Eq. (la) separately for the absolute value \tp\ and the 
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phase (j) of the order parameter 



dt 

'dt 



1^1(1 -IV'I' 

1 djn 



uj-i/ip ds 



(4a) 
(4b) 



Here s is the arc-coordinate along the ring and we used 
the condition div(js + jn) = 0. It is obvious from Eq. 
(4a) that if the RHS of Eq. (4a) is negative jV'l decreases 
in time and if the RHS is positive {ipl wih increase in 
time. Because the second derivative of \ip\ is always pos- 
itive (at least near the phase-slip center - see Fig. 4(a)) 
the different time dependence of lipl is governed by the 
term — p^lV'l- From Fig. 4(b) it is clear that after the 
phase slip the value of p is less than before this moment, 
with practically the same distribution of 1-01 • It is this 
fact which is responsible for an increase of the order pa- 
rameter just after the moment of the phase slip. But at 
some moment of time the momentum p can become suf- 
ficiently large, making the RHS of Eq. (4a) negative and 
as a consequence starts to decrease. 
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FIG. 4: Distribution of the absolute value of the order param- 
eter (a) and gauge-invariant momentum (b) near the phase 
slip center at different moments of time for a ring of i? = 10^ 
and with u = 3. The phase slip occurs at t ~ 2.7tgl- 

Based on our numerical calculations we may state that 
for every value of the order parameter in the minimal 
point, there exists a critical value of the momentum p™"' 
such that if the value of the momentum in this point is 
less than the order parameter will increase in time. 
In the opposite case the order parameter decreases which 
leads ultimately to the phase slip process. Therefore, 
when after a phase slip p increases fast enough such that 
at some moment the condition > p™*" is fulfilled, 

the order parameter will start to decrease which leads 
than to a new phase slip process. 

From our results we can conclude that the change of 
p (or phase of the order parameter) with time and of 
with time has a different dependence on u. Indeed 
from Figs. 2 and 3 it follows that with increasing u the 



time relaxation of IV'I becomes larger with respect to the 
relaxation time for p (for example, the amplitude of os- 
cillations of is decreased). Moreover, the relaxation 
time for p depends not only on u but also on the history 
of the system: the larger the number of phase slips which 
have occurred in the system the longer the time becomes 
between the next two subsequent phase slips. In order to 
understand such behavior we turn to Eq. (4b). 

The numerical analysis shows that the second term in 
the RHS of Eq. (4b) is only important for some length 
near the phase slip center. This length is nothing else 
than the length over which the electric field penetrates 
the superconductor. We checked this directly by solving 
Eqs. (la,b) for this situation. But this length is the decay 
length \q of the chargpJmbalance in the superconductor 
(see for example Refs.E^Lj). We numerically found that 
\q varies with u (in the range u = 1 — 100) as Xq ^ 
y-o.27_ means that the relaxation time of the charge 
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imbalance is tq 

Lets take the derivative d/ds on both side of Eq. (4b) 
and integrate it over a distance of Xq near the PS center 
(far from it we have jn 0) we obtain the Josephson 
relation (in dimcnsionlcss units) 



dA<j> 
dt 



(5) 



where Acj) = </)(-!- Aq) — 0(— Ag) is the phase difference 
over the phase slip center which leads to the voltage 
V — —{{p{+Xq) — ip{~Xq)) and where jn(0) is the nor- 
mal current (or electric field in our units) in the point of 
the phase slip. From Eq. (5) it follows immediately that 
the relaxation time for the phase of the order parameter 
near the phase slip center is 



1 



T4, 



Aq < |jn(0)| >' 



(6) 



where < • > means averaging over time between two con- 
secutive phase slips. This result allows us to qualitatively 
explain our numerical results. Indeed, from Fig. 2(c) it 
is apparent that < |jn(0)| > decreases after each PS and 
as a result it leads to an increase of r^. With increas- 
ing u, Xq decreases and increases. Fitting our data 
(see Fig. 5) leads to the dependencies ~ and 
T0 ~ u^ '^^ for rings with radius i? = 10 and R = 15, 
respectively (here r<ji was defined as the time between 
the first and the second phase slip). This is close to the 
expected dependencies which follow from Eq. (6) and 
from the above dependence of Xq(u) (quantitative dif- 
ferences follows from uncertainty in our finding and 
Xq). Besides, Eq. (6) allows us to explain the decrease 
of with increasing radius of the ring (see inset in Fig. 
5). Namely, during the time between two phase slips the 
gauge-invariant momentum decreases as Ap ~ 1/R (be- 
cause V(/) increases as ~ 27r/27ri?) in the system. Far 
from the phase slip center the total current practically is 
equal to js ~ P- Because div(js + jn) = in the ring we 
directly obtain that during this time the normal current 
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density in the point of the phase shp also decreases as 
~ 1/R. Taking into account Eq. (6) we can conclude 
that should vary as ~ 1 /R (at least for a large radius 
and for the first phase slip). The behavior shown in the 
inset of Fig. 5 is very close to such a dependence (it is 
interesting to note that in contrast to the time 
does practically not depend on R). 

On the basis of our results we can make the following 
statement: when the period of oscillation (time of change 
of the phase of the order parameter) becomes of the or- 
der, or larger, than the relaxation time of the absolute 
value of the order parameter the next phase slip becomes 
impossible in the system. 

This result can be applied to the system of Ref.EZi where 
a current carrying wire was studied. The authors found 
a critical current jc — 0.335 for u = 5.79 and jc = 0.291 
for u — \2. From our observation follows that the phase 
slip solution may be realized as a stable solution when 
T|v,| > ^ 1/(Aq < |jn(0)| >). We found that Ag 

T|^l ^ u"-^ and < |jn(0)| >~ jext which leads 
to the critical current jc ~ u^^'^^ which decreases with 
increasing 

Above we found that the time scale governing the 
change in the phase does not coincide with the relax- 
ation time of the absolute value of the order parameter. 
This difference is essentially connected with the presence 
of an electrostatic potential in the system. In order to 
demonstrate this we performed the following numerical 
experiment. We neglected in Eq. (la) and found that 
the number of PS is now independent of u and R, and AL 
equals unity. Moreover, it turned out that the time scale 
to is an order of magnitude larger. This clearly shows 
that the electrostatic potential is responsible for the ap- 
pearance of a second characteristic time which results in 
the above mentioned effects. 

In an earlier papeilij the question of the selection of 
the metastable state was already discussed for the case 
of a superconducting ring. However, these authors ne- 
glected the effect of the electrostatic potential and found 
that transitions with AL > 1 can occur only when the 
magnetic field— (called the induced electro-moving force 
(emf) in Ref.Ej) increases very quickly. In their case 
these transitions were connected to the appearance of 
several nodes in if) along the perimeter of the ring and 
in each node a single PS occurred. We reproduced those 
results and found such transitions also for larger values 
of Ai?. However, simple estimates show that in order 
to realize such a situation in practice it is necessary to 
have an extremely large ramp of the magnetic field. For 
example, for Al mesoscopic samples with ^(0) = lOOnm 
and i? = 10^ the corresponding ramp should be about 
10'^ — lO'^T/s. With such a ramp the induced normal 
currents in the ring are so large that heating effects will 
suppress superconductivity. 

The transitions between the metastable states in the 
ring are not only determined by u and i?, but also, the 
presence of defects play a crucial role. Their existence 
leads to a decrease of Ai. This is mainly connected 
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FIG. 5: Dependence of the initial nucleation time ta (curves 
1,2), the relaxation time of the order parameter r|^|, (which 
for definiteness we defined as the time of variation of the order 
parameter from 0.4 till the first PS - curves 3,4) and the relax- 
ation time of the phase of the order parameter (we defined 
it here as the time between the first and second PS - curves 
5,6) are shown as function of the parameter u. Dotted(sofid) 
curves are for a ring with radius R = 10(15). Dashed curve 
corresponds to the relaxation time of the charge imbalance 
(as obtained from the expression tq — 5.79Aq). In the inset, 
the dependency of on the ring radius is shown for different 
values of the parameter u. 



to the fact that with decreasing pc (as a result of the 
presence of a defect) the value of < |in(0)| > decreases 
even at the moment of the first phase slip (because < 
bn(0)| > cannot be larger than < |js(0)| >— Pc) and 
hence increases. We can also explain this in terms of 
a decreasing degeneracy of the system. For example, if 
we decrease Pc with a factor of two (e.g. by the presence 
of defects) it leads to a twice smaller value for <I>c and 
Leq. But our numerical analysis showed that the effect 
of defects is not only restricted to a decrease of Pc and 
Leg. To show this we simulated a defect by inserting an 
inclusion of another material with less or zero T^- This is 
done by inserting in the RHS of Eq. (la) an additional 
term p{r)'ip where p{r) is zero except inside the defect 
where p{r) = a < 0. The magnetic field was increased up 
to He of the ideal ring case. We found that the number of 
PS was smaller as compared to the case of a ring without 
defects. The calculations was done for defects such that 
Pc was decreased by less than a factor of 2 in comparison 
to the ideal ring case. In contrast, a similar calculation 
for a ring with nonuniform thickness/width showed that, 
even for 'weak' nonuniformity (which decreases Pc by less 
than 20%), AL was larger than for the ideal ring case 
and the final vorticity approaches Leg. This remarkable 
difference between the situation for a defect and the case 
of a nonuniformity may be traced back to the difference 
in the distribution of the order parameter: even in the 
absence of an external magnetic field the defect leads to 
a nonuniform distribution of the order parameter which 
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is not so for the nonuniform ring case. A more thorough 
study of the effect of defects will be presented elsewhere. 



IV. SUPERCONDUCTING RING WITH 
NONZERO WIDTH 

All the above results were based on a one-dimensional 
model which contains the essential physics of the decay 
and recovery of the superconducting state in a ring from 
a metastable state to its final state. However, even in 
the case when the effect of the self-field can be neglected 
the finite width of the ring may still lead to important 
additional effects (for example, a finite critical magnetic 
field). In order to include the finite width of the ring 
into our calculation we considered the following model. 
We took a ring of mean radius R — 12^, width 3.5^, and 
thickness less than ^ and A. These parameters arei-close 
to those of the ring studied experimentally in Ref.a and 
they are such that we can still neglect the self-field of the 
ring. The obtained magnetization curves of such a ring 
are shown in Fig. 6 as function of the magnetic field and 
for two values of u. Those results were obtained from 
a numerical solution of the two-dimensional Ginzburg- 
Landau equations, Eqs. (la,b). The magnetic field was 
changed with steps /S.H over a time interval which is 
larger than the initial part of the phase slip process 
(for our parameters this procedure leads practically to 
an adiabatic change of the magnetic field) . 

From Fig. 6 we notice that the value of the vorticity 
jumps, AL, depends sensitively on the parameter zJla. 
We found that the phase slips occur in one particular 
place along the perimeter of the ring. However, in con- 
trast to our previous one-dimensional case, AL depends 
also on the applied magnetic field. The reason is that 
for a finite width ring the number of metastable|-s±ates 
decreases with increasing magnetic field (see Ref.cj). It 
means that the system cannot be moved far from equilib- 
rium with a large superconducting current density (be- 
cause the order parameter is strongly suppressed by the 
external field) at high magnetic field. Hence, the value of 
< |jn(0)| > will be much smaller in comparison with the 
one at low magnetic fields and is larger or comparable 
with r|^| even for high values of u and even for the first 
PS. Thus, the effect of a large magnetic field in case of a 
finite width ring is similar, to some respect, to the effect 
of defects for a one-dimensional ring. 

Our numerical results are in qualitative agreement with 
the expprimcntal results of Pedersen et al. (see Fig. 2 
of Ref.cl). Unfortunately, no quantitative comparison is 
possible because of a number of unknowns, e.g. the value 
of pc (it is necessary to have the dependence of M{H) as 
obtained starting from zero magnetic field), the value of 
^ is not accurately known and hence the ratio i?/^ can 
therefore only be estimated. The value of both these 
parameters have a strong influence on the value of the 
vorticity jumps AL. 

Let us discuss now finite temperature and screening 
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FIG. 6: Dependence of the magnetization of a finite width 
ring on the external magnetic field for two values of the pa- 
rameter u. The results are shown for a sweep up and sweep 
down of the magnetic field. 

effects (i.e. finite A). Nonzero temperature leads to the 
possibility of a perturbative appearance of phase slip cen- 
ters at U* < He- But if the field H* does not differ too 
much from He the average < |in(0)| > is not altered 
strongly. As a result, after the appearance of the first 
thermo-activated PS center the second phase slip may ap- 
pear automatically (if the ratio T|0|/t0 is large enough) 
which leads to an avalanche- type of process. Unfoctu- 
nately, it is impossible to apply the results of RefBd in 
order to calculate the probability of this, process even for 
the first PS. The reason is that in Ref.Eil the effect of the 
electrostatic potential (which play a crucial role as was 
mentioned above) was neglected. 

A finite A (and large enough width of the ring) leads to 
a nonuniform distribution of the momentum p over the 
width of the ring which can play an important role on the 
above considered effects. As we mentioned above such a 
finite width ring may be modelled as a strip with Itans- 
port current in an external magnetic field. In RefEj the 
condition for the entry of the first vortices in a narrow 
superconducting strip was studied. It turned out, that 
the period of vortex chain entry depends essentially on 
the distribution of p over the width of the strip. When 
this distribution is uniform the period is infinite. But 
even small nonuniformities in this distribution leads to 
a finite period. If we translate their results to our ring 
geometry, the period of the vortex chain must now be dis- 
crete 2Tr R/n (n = 1, 2, ...). Under a certain condition the 
entry of vortices becomes possible, not through a single 
point along the perimeter, but through two, three and 
more points (without increasing the ramp of the mag- 
netic field). A finite A only increases this effect because 
it leads to larger nonuniformities in the dependence of p 
over the ring width. In this situation there is a competi- 
tion between the process of the appearance of additional 
nodes along the ring perimeter and the number of PS in 
these points. At high magnetic fields these effects be- 



7 



comes negligible small (because superconductivity only 
exists near the edges and the ring effectively can be con- 
sidered as two one-dimensional rings) and we will have 
the situation as discussed in the present work, and, re- 
cently incj. 

V. CONCLUSION 

In conclusion, we studied how an unstable supercon- 
ducting state of a superconducting ring evolves in time 
and transit to its final state. The latter is not neces- 
sary the thermodynamic equilibrium state and may be 
another metastable state with a different vorticity but 
which is stable in time. The transition between the dif- 
ferent superconducting states occurs through a phase slip 
center which is a point along the ring where the supercon- 
ducting amplitude decreases to zero abruptly resulting in 
the change of the vorticity of the superconducting state 
with one unit. The waiting time, or the creation time for 
the first phase slip, is found to be two orders of magni- 
tude larger than the subsequent time intervals between 



consecutive phase slips. The latter time is connected with 
the time relaxation of the charge imbalance in the super- 
conductor and increases the closer the system becomes 
to the final state. This circumstance allows to find this 
time (also as the relaxation time of the absolute value 
of the order parameter) from magnetic measurements of 
superconducting rings of large radius. Our theoretical 
findings are in agreernent with recent experimental re- 
sults of Pedersen et aln. 
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